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ABSTRACT 



Advances in hydrodynamical simulations have pro- 
vided new insights into the effects of convection on 



the frequencies of solar oscillations. As more accurate 
observations become available, this may lead to an 
improved understanding of the dynamics of convec- 
tion an d the interaction betw een convection and pul- 
sation ( Rosenthal et al. 1999] ) . Recent high-resolution 
observations from the SUl/MDI instrument on the 
SOHO spacecraft have provided the so-far most- 
detailed observations of high-degree modes of solar 
oscillations, which are particularly sensitive to the 
near-surface properties of the Sun. Here we present 
preliminary results of a comparison between these 
observations and frequencies computed for models 
based on realistic simulations of near-surface con- 
vection. Such comparisons may be expected to help 
in identifying the causes of the remaining differences 
between the observed frequencies and those of solar 
models. 



including the role of correlations, which result in a 
mean structure which does not satisfy the locally- 
valid equation of state, and radiative effects due to 
opacity inhomogeneities, which result in hotte r sur 



face layers for a giv en effective temperature (see Stein 



fc Nordlund 1998a) 



Helioseismology has considerable potential as a tool 
for testing the accuracy of numerical models of these 
surface regions through comparison of predicted and 
measured frequencies. A major contribution is the 
effect of turbulent pressure on the mode frequencies 
( Kosovichev 1995] ). In such a comparison, it is vi- 
tal to distinguish between model effects on frequen- 
cies, resulting from errors or inaccuracies in the mean 
structure, and modal effects resulting from uncer- 
tainties in the physics of wave propagation in these 
turbulent layers. We use numerical simulations of the 
surface layers to handle model effects while studying 
various physically-motivated assumptions about the 
modal physics. 



Key words: solar oscillations; convection. 



2. DATA ANALYSIS 



1. INTRODUCTION 



In stellar-structure theory, the standard Mixing- 
Length Model of convection (MLT) continues to be 
the industry standard, despite its many limitations. 
For example MLT cannot account for the broaden- 
ing of photospheric spectral lines, presumably due to 
overshooting, nor can it model the effect of turbulent 
pressure on the mean structure in a consistent way. 
More subtly perhaps, the reduction of the true 3-D 
radiation hydrodynamics to 1-D hydrostatic stratifi- 
cation eliminates some significant contributing effects 



Our data analysis procedure is described in more de- 
tail by Rhodes et al. (1997, 1998). A 60.75 day 
time series of MDI Full-Disk Dopplergrams, begin- 
ning May 23, 1996, was used. After limited gap- 
filling, the final duty cycle was 97.3 %. The so- 
called "averaged-spectrum" method is used to ob- 
tain frequency estimates, with results from low- and 
intermediate-degree power spectral peaks being used 
to obtain a correction due to asymmetry which is 
then applied to the frequencies obtained by ridge- 
fitting at higher degrees. Results are obtained for 
degrees up to t — 1000. 
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3. THE STANDARD SOLAR MODEL 



The standard solar model we use is 



Christenscn-Dalsgaard et al. (1996) 



Model S of 
It uses the 



UJr'AL equation ot state ( Rogers ct al. 1996), OPAL 



opac ities ( Iglcsias ct al. 1992] ) at high temperature 
and (Kurucz~1991 ) opacities m the atmosphere. The 
atmosphere is determined by a T — r relation (where 
T is temperature and r is optical depth) fitted to the 



Our philosophy is to construct numerical simulations 
and standard MLT envelope models each using the 
same physics, thus allowing us to isolate the effects 
of 3-D convective dynamics. In this case we used 
simul ations calculated with th e MHD equation of 
state ([Hummer fc Mihalas 1988|; |Mihalas ct al. 1988 



HSRA reference at mosphere (Gingcrich ct al. 1971 
and standard MLT ( Bohm-Vitcnsc 195£ ) is used for 
the convection zone. The effects of turbulent pressure 
are not included in Model S. 



Figure 1 shows a comparison between adiabatic 
eigenfrequencies of Model S and the data. The fre- 
quency residuals have been scaled with the quantities 
Q n i which represent the ratio of the mode inertia of 
a particular mode to that of a radial mode of the 
same frequency. Mode inertias were normalised at 
the photosphere, defined as the location where the 
temperature is equal to the effective temperature. It 
is evident that, except for modes of low order, the 
scaled residuals are essentially a function of frequency 
alone, an indicator that the predominant effect re- 
sponsible for the residuals lies close to the solar sur- 



Dappen et al. 19"8£~ Mihalas et al. 199C| ) and Kurucz 
opacities. The effective temperature of the simula- 
tion was 5777K. 

A standard envelope model was constructed with the 
same physics, and using standard MLT with the ad- 
dition of a parameterised turbulent pressure of the 
form pt — Ppv 2 , where p is density and v is the the 
convective velocity, obtained from the mixing-length 
formulation. The extra parameter (3 and the usual 
mixing length parameter were then adjusted to al- 
low a continuous matching of density, pressure and 
temperature. Because the matching point is at a 
depth where the turbulent pressure is already very 
small (of order 1% of the total pressure) we do not 
consider this matching process to be a useful pa- 
rameterisation of the relation between mixing-length 
velocity and turbulent pressure in general. Strik- 
ingly, the resulting matched model has a convec- 
tion zone depth of 0.286-R©, very clos e to the he 



fade~((Jhristcnsen-Ualsgaarcl & Berthomieu 1991J 



lioscismic value of (0 . 287 ± 0.003) .R (|Christcnscn 
Dalsgaard et al. 1991; Kosovichcv k FccTorova 1991)' 



or (0.287 ±O.UO1)^ (Easy fc Antia 1997). This is a 
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Figure 1. Measured frequency residuals in the sense (ob- 
servations) - (model), scaled by Q n i for selected modes 
in the range < £ < 1000. The computed frequencies 
are for Model S of Christens en- Dalsgaard et al. (1996). 
Radial orders of the lowest-order modes are indicated f, 
pi, andp2- 



MODEL CONTRIBUTIONS AND 
NUMERICAL SIMULATION 



remarkable, parameter-tree validation ot the simula- 
tions. 



5. MODAL CONTRIBUTIONS 



The process described above yields a hydrostatic 
pressure and density stratification for a patched 
model based on the simulation. We consider only 
adiabatic oscillations of this patched model. Under 
the assumption of adiabaticity, the modal effects are 
essentially confined to the relation between the La- 
grangian perturbations Sp and Sp in pressure and 
density, 



5p 
P 



Sp 
P 



(1) 



When convection is ignored, T x = (d\np/d\n p) s , 
the derivative being at constant specific entropy s, is 
simply the thermodynamic adiabatic exponent. Here 
we incorporate modal effects by the use of different 
recipes for Ti; in particular, they define the adiabatic 
sound speed c = (Tip/p) 1 / 2 . We consider two pos- 
sibilities: the Gas Gamma Model (GGM) and the 
Reduced Gamma Model (RGM). 

In the GGM we simply use the horizontal and tem- 
poral mean value of the thermodynamic T% from the 
simulations. In the RGM, this is replaced by the so- 
called "reduced gamma" : 



We treat the model contributions to the residuals by 
using radiation-hydrodynamic simulations of the sur- 
face layers to construct new mean models, by patch- 
ing temporal and spatial averages of the simulations 
to convective envelope models constructed with stan- 
dard physics. The simulations are described in detail 
by Stein & Nordlund (1989, 1998). 
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where p g is the gas pressure and p is the total pres- 
sur e. The reduced gamma was originally introduced 
by Rosenthal et al. (1995), who argued that if the 
Lagrangian perturbation to the turbulent pressure, 



3 



5pt, were zero then 
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(3) 



Rosenthal et al. (1999) show that the definition (|2[) 
is the correct one to use if Sp t — provided certain 
other conditions are also met, principally that the 
time variation of the convective energy fluxes can be 
neglected. 

We also construct a Standard Envelope Model (SEM) 
which uses the same physics as the simulations, has 
no turbulent pressure, and is calibrated to have the 
same convection-zone depth as the patched models. 
The SEM is constructed using the Rosseland Mean 
of the opacities used in the simulations and a T — t 
relation determined from the averaged simulation it- 
self. Thus the differences between the models should 
isolate truly convective effects. 



6. MODEL DIFFERENCES 



Christcnsen-Dalsgaard & Thompson (1997) have 
shown that when model differences are expressed in 
Lagrangian terms (i.e., on a fixed mass scale) the fre- 
quency change due to near-surface perturbations can 
be written approximately as 



Su>ri 



KJr) dr 



(4) 



where S m v/v is the Lagrangian perturbation to the 
quantity v = Ti/c, and where the kernels K™ l c (r) 
can be calculated from the structure and eigenfunc- 
tions of the reference model. For an isothermal layer, 
v = 2u c /g, where u c is the acoustic cutoff frequency. 
Hence changes in v correspond physically to changes 
in the size of the acoustic cavity. 

Figure 2 shows Lagrangian differences in v between 
the various models. The GGM-SEM differences 
largely reflect changes in density due to the elevation 
of the photospheric regions by turbulent pressure. In 
addition, the RGM-SEM differences show the effect 
of the turbulent reduction of Ti. Comparing (GGM- 
SEM) with (GGM-Model S) we can see that the dif- 
ferent physics between SEM and Model S (principally 
OPAL v. MHD equation of state and a different T—t 
relation) produces further small changes confined ex- 
tremely close to the surface. 



7. OSCILLATION FREQUENCIES 



We restrict ourselves to modes trapped within the 
convection zone. This allows us to include direct 
comparisons also with Model S. Figure 3 shows that 
the frequency differences between GGM and SEM are 
similar in magnitude and shape to those between the 
data and Model S (Figure 1), whereas corresponding 
differences between RGM and SEM are much larger 
than required. GGM-Model S differences are broadly 
similar to GGM-SEM differences. 
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Figure 2. Logarithmic Lagrangian ( at fixed mass ) differ- 
ences inv = Y\j c between patched and comparison enve- 
lope models (solid and dashed lines): SEM is the compar- 
ison model, using mixing-length theory, while GGM and 
RGM are patched models assuming the averaged thermo- 
dynamic Ti and the reduced T\, respectively. Also shown 
is the difference between the GGM envelope and standard 
solar model (Model S). 



Figure 4 is a direct comparison between GGM fre- 
quencies and the data. The residuals are much re- 
duced compared to Figure 1. However the "bump" 
below 3mHz suggests that the combined model and 
modal differences in Figure 2 are too deeply penetrat- 
ing. Thus the Gas Gamma hypothesis cannot be a 
completely accurate description of what is occurring. 



INTERPRETATION AND CONCLUSIONS 



Our statement that the remaining discrepancy in fre- 
quency must be due to modal, rather than model, 
effects is based on the belief that the simulations are 
fundamentally correct. Why do we have such faith 
in the simulations? The principal model effect on 
the frequencies is the elevation of the photospheric 
layers by the combined effect of turbulent pressure 
and an effect of the 3-D dynamics, whereby hotter 
regions, which have higher opacity contribute less to 
the emitted radiation. The higher mean temperature 
of the 3-D model is thus "hidden" from view, but is 
reflected in the pressure stratification. 

We have carried out a number of convergence tests on 
the simulations involving runs at various resolutions. 
The thermal stratificat ion is rather insensitive to nu- 
merical resolution (cf. Stein & Nordlund 1998b, Fig. 
26). The variation of turbulent pressure with resolu- 
tion is shown in Figure 5. 



The highest-resolution runs already produce the cor- 
rect turbulent broadening, asymmetries, and blue- 
shift of spectral lines and so the maximum turbulent 
pressure cannot i ncrease much beyond the largest val- 



ues found so far ( Nordlund & Dravins 1990 : pravins 
k Nordlund 1990a[rpravins k NoraTunrtH$90E|j7Tri 



tact, despite the apparent steady increase of turbu 
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Figure 3. Scaled differences between frequencies of 
patched (RGM and GGM) and comparison (SEM) enve- 
lope models and between the GGM envelope and the stan- 
dard solar model (Model S). 
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Figure 5. Mean turbulent pressure divided by total (gas 
+ turbulent) pressure as a function of depth for four dif- 
ferent runs, with different numerical resolutions: 253 2 x 
163, 125 2 x 82, 63 2 x 63, and 32 2 x 41. (Note that 
the total duration of the 253 2 X 163 run is only 10 min- 
utes, while we estimate that a minimum duration of 30-60 
minutes is required for statistically stable results.) 
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Figure 4- Measured frequency residuals in the sense (ob- 
servations) - (model), scaled by Q n l- The model frequen- 
cies are for the GGM, with the (unreduced) gas Ti. 



lent pressure with numerical resolution shown in Fig- 
ure 5, the density structure varies very little between 
the various runs. Hence, we expect that the elevation 
of the atmosphere is nearly fully accounted for at our 
highest resolutions. 

What then are the nature of the modal effects which 
account for the rema ining frequency differences? One 
clue is provided by Btein & Nordlund (1991), who 
show that the time variation ot the turbulent pres- 
sure is related to the compression in a much more 



complicated fashion than that implied by either the 
GGM or the RGM. In practice the turbulent compo- 
nent of gamma, 91npt /<91np, can be negative, and is 
a strong function of height. Evidently, we need to use 
such information from the simulations themselves to 



determine a better model for Ti than the two simple 
cases considered here. 

However, the limitations even of this approach are 
evident from consideration of the f modes. These are 
almost entirely insensitive to the mean solar struc- 
ture and so no treatment based on variations in either 
Ti or the hydrostatic structure will ever improve the 
agreement between measurement and theory in their 
case. The principal contribution to the f-mode resid- 



bv turbulent motions (Murawski & Roberts 1993a 


Murawski & Roberts 1993b 




(Jruzmov 499i J, an ct- 



realistically. However, Murawski et al. (1998) have 
derived a dispersion relation for the f mode in a sim- 
ple model containing a random velocity field, and 
shown that the frequency shift and linewidth of the f 
mode for this model are consistent with the proper- 
ties obtained from the high-resolution MDI data. It 
is likely that turbulent advection results in frequency 
shift of the p modes in addition to the turbulent pres- 
sure effect. It may be, once again, that we will have 
to turn to the simulations to provide us with clues as 
to how to describe these processes. 
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